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Abstract 

■ For the statistical multifragmentation model the critical indices a , /3, 7 , are cal- 

culated as functions of the Fisher parameter r. It is found that these indices have 
different values than in Fisher's droplet model. Some peculiarities of the scaling rela- 
tions are discussed. The basic model predicts for the index r a narrow range of values, 
1.799 < t < 1.846, which is consistent with two experiments on nuclear multifragmen- 
tation. 
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r-| '. I. Introduction 



The statistical multifragmentation model (SMM) Jl], |2|] rather successfully explains the 
multiple production of intermediate mass fragments |3], |J. On the other hand it is also 
an example of a statistical model exhibiting a phase transition (PT) of the liquid-gas type 
|4], [J. Therefore, the SMM was studied extensively in order to clarify not only the relationship 
between multifragmentation and the nuclear liquid-gas PT, but also to elucidate the connection 
of the SMM to both Fisher's droplet model (FDM) @ and the lattice gas model [0, [§. For 
the last purpose a simplified version of the SMM was suggested || [3~Qf| . 

In recent works pi |12| an exact analytical solution of a simplified version of the SMM 
was found in the thermodynamic limit. This model is mathematically very similar to the 
original FDM. However, it differs from the FDM in two aspects: (i) the nonzero volume of 
the nucleons is taken into account, and (ii) the free surface energy is non- negative. These 
two main differences generate distinct results concerning the properties of the critical point: 
the SMM predicts the existence of a critical point for the parameter r < 1 and a tricritical 
point for 1 < r < 2 (see below), whereas for the FDM this parameter is assumed to be 
larger, r > 2. Such a difference is not surprising since the FDM considers the droplets to have 
zero volume, i.e., its validity is restricted to low particle densities. Consequently, the FDM 
conclusions regarding the critical behavior must be considered as mainly suggestive in relation 
to the behavior of real fluid systems or of more realistic models. Therefore, it is necessary to 
study how the critical indices of the SMM relate to those of the FDM. It will be shown that 
the critical indices in the SMM and the FDM have different values, i.e., these models belong 
to different universality classes. 

Another important issue is the comparison of the critical indices of the SMM with the 



corresponding values found in numerical studies |L3[ and in experiments. As we shall show 



the experimental data for the index r are consistent with our predictions for the SMM. 



The paper is organized as follows. In Sect. II a brief description of the simplified version 
of the SMM is given. In Sect. Ill some useful formulae are derived and the critical indices of 
this model are found. The scaling laws and the relations between the critical indices of the 
SMM are discussed in Sect. IV. The conclusions are given in Sect. V. 

II. The Statistical Multifragmentation Model 

~1| a simplified version of the SMM is solved analytically with the method 
The interaction of the fragments is taken into account by replacing the total 



In Refs. [11 



suggested in JL4 



volume V by the free volume Vf = V — bA, where A is the number of nucleons in the system 
and b is the eigenvolume of a nucleon. In the grand canonical ensemble the pressure of the 
system is found for the liquid and gaseous phase 

li + W{T) 
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p g (T,fi) = TF(p g ,T, fi) 



respectively, where 



F(p g ,T,n) = g{T) 
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and 



with g(T) = [mT/(2vr)] 3/2 . The mass of a nucleon is denoted as m, z\ = 4 is the number 
of spin-isospin states of a nucleon, v = p + W(T) — bp g is the shifted chemical potential 
of the gaseous phase and W(T) is the free energy per particle inside of a fragment. The 
critical temperature of the system is denoted as T c and e = (T c — T)/T c defines the reduced 
temperature. The actual parameterization of W(T) and the free surface tension a(e) can be 
found in Ref. [11, 12]. In this paper it is assumed that W(T) and all its derivatives are regular 
functions of temperature and that the free surface tension obeys a(e) = a e^ for small positive 
values of e and vanishes for e < 0. The surface area of a fragment of k nucleons is taken in the 
general form k a . In the SMM one usually chooses £ = 5/4 and o = 2/3, but in the following 
treatment £ > 1 and < a < 1 will be regarded as free parameters. 

The particle densities in the liquid and in the gaseous phase are p\ = (dpi/dp) T = 1/6 and 
p g = {dp g /dp) T = pa/ (1 + bpid), respectively. The density 
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(3) 



Pid(T,n) = g(T) 



v-W(T) 



+ Ei(£, v) 
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has the meaning of the ideal gas density in the limit 6^0. 

When both pressures coincide, pi = p g , v = , the liquid-gas PT occurs. This equality 
defines the PT equilibrium curve p*{T) in the (T, /x)-plane. For T < T c the surface tension 
a(e) is nonzero and the PT is of 1st order, since p id (T,p) < oo and consequently p g < p% on 
the curve p*(T). The critical point at T = T c corresponds to the critical density p c — p t — 1/b. 
Since at the critical point v — and the surface tension vanishes, it follows from Eqs. (|l|) and 
(0) that for r < 1 the function T is divergent and, therefore, p c is infinite. For 1 < r < 2 one 
finds similarly that p c < oo , p c = 1/6 and a 2nd order PT for T > T c (see Ref. []i~l , |12| for more 
details). For r > 2 on the phase equilibrium line both ^(^(T), T, p*(T)) and pa(T, pi*(T)) 
are finite for any finite temperature, i.e., the mixed phase region extends formally to infinite 
temperatures. Therefore, this case exhibits no critical point at all in contrast to the FDM and 
is not considered in the following. 



III. The Critical Indices of the SMM 

The critical exponents a', j3 and 7' describe the temperature dependence of the system 
near the critical point on the coexistence curve v — 

/ \e\~ a , fore<0, 
C > ~ \e-«, for £ >0, (5) 

Ap ~ , for e > , (6) 

k t ~ e" 7 ' , for £ > , (7) 

where Ap = pi — p g defines the order parameter, c p = ^ (^f) denotes the specific heat at 

fixed particle density and kt = 1 (ff) T * s ^ ne isothermal compressibility. The shape of the 
critical isotherm for p < p c is given by the critical index S 

Pc-P ~ (Pc - p) 5 for e = . (8) 

The tilde indicates e = hereafter. 

The calculation of a and a' requires the specific heat c p . With the formula [15 



Cp (T,p) \(d 2 p\ (d 2 p\ 



T p \dT 2 J \dT 2 



p 



one obtains the specific heat on the PT curve by replacing the partial derivatives by the total 
ones [ff6| . The latter can be done for every state inside or on the boundary of the mixed phase 
region. For the chemical potential p*(T) = bp*(T) — W(T) one gets 

1MH^ + ™ (io) 

where the asterisk indicates the condensation line (v = 0) hereafter. Fixing p = p c = pi = 1/b 
one finds c* ; (T) = T d and, hence, obtains similarly for both a and a' 

a = a' = . (11) 

To calculate (3 and 7' the behavior of the series 

00 

£ q M) = k9 ~ T z- M ' k ° (12) 

k=2 

for small positive values of e should be analyzed (A = a Q /T c ). In the limit e — > the function 
E q (e, 0) remains finite, if r > g + 1, and diverges otherwise. For r = q + 1 this divergence is 
logarithmic. The case r < q + 1 is analyzed in the following. 

With the substitution z^ = k (As^) ^ it follows 

r-q 00 

S q (s,0) = (^) CT £ 4" T e - ^ ■ (13) 

k=2 



Since Azk = ( y 4 £< ') ^ = the series on the r. h. s. of fll3|) converges to an integral for e — > 



T-q-l -- ct 

q (e,0) = (Ae c ) CT / ds 2«- T . 



(14) 



The assumption q — r > — 1 is sufficient to guarantee the convergence of the integral at its 
lower limit. Because of the finite value of the above integral one obtains the following result 
for e — ► 

i(r-?-i; 



S q (e,0)~ i n | £ | ; 
For the calculation of 5 one obtains similarly 



, if r < q + 1 , 
if r = g + 1 , 
if r > g + 1 . 



(15) 
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S q (0,z>) ~ { ln|z>|, 



if r = q + 1 , 
if t > g + 1 . 



(16) 



From Ap = ~ 1+ l p . d h follows Ap ~ p* d 1 for e — > . Therefore, inserting q = 1 into Eq. 
]) immediately gives /3 



/? = £(2-r) 



(7 



From p g = Pid/(1 + fepid) an d Eq. (|4]) it follows that 



"-wen ^ , N 
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(17) 
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Using Eq. (|i"5"D one easily finds 7' from Eq. (|18|) 



y = - 

cr 



r "2 



(19) 



Near the critical point it follows p c — p ~ l/pid- With (||) one finds using Eq. ([TBI) Pc — P 
[£i(0, vyf 1 ~ i> 2 ~ r . Substituting 1 = e^ k e^ k yields for small values of v 



Pc-P = g(T c )T c 



00 
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Combining the result above with the expression found for p c — p one obtains the critical index 
5 

T - 1 

8=1 < (21) 



2 - r 



which is independent of the choice of ( , o and r. 



IV. Scaling Relations of the SMM 

In the special case £ = 2a the well-known exponent inequalities proven for real gases by 



Fisher [17] 
Griffithsra 



Liberman|L9| 



a' + 2(3 + i > 2 , (22) 
a' + /3(l + 8) > 2 , (23) 
i + P(l-5) > , (24) 



are fulfilled exactly for any r. (The corresponding exponent inequalities for magnetic systems 
are often called Rushbrooke's, Griffiths' and Widom's inequalities, respectively.) For £ > 2a, 
Fisher's and Griffiths' exponent inequalities are fulfilled as inequalities and for £ < 2a they 
are not fulfilled. The contradiction to Fisher's and Griffiths' exponent inequalities in this last 
case is not surprising. This is due to the fact that in the present version of the SMM the 
critical isochore p = p c = p\ lies on the boundary of the mixed phase to the liquid. Therefore, 
in expression (2.13) in Ref. |T7| for the specific heat only the liquid phase contributes and 
the proof of Ref. [[17]] following (2.13) cannot be applied for the SMM. Thus, the exponent 
inequalities fl2|) and fl23|) have to be modified for the SMM. Using Eqs. (JT^), (0), © and 
(pl|) one finds the following scaling relations 

a' + 2(3 + i = - and a' + (3(1 + 6) = - . (25) 
a a 



Liberman's exponent inequality (|24|) is fulfilled exactly for any choice of C and a. 

Since the coexistence curve of the SMM is not symmetric with respect to p — p c , it is 
interesting with regard to the specific heat to consider the difference Ac p (T) = c* (T) — c* ; (T), 
following the suggestion of Ref. JRJ. With Eq. ([XQ|) and noting that 1/p* —b = 1/ p* d it follows 



A ^ T) = i£m-dfr- ■ (26) 

The most divergent term in Eq. (^) yields for ( > 1 

f 2 - £ , ifr<a + l, 



Then it is a' s > for Q/a < 2. Thus, approaching the critical point along any isochore within 
the mixed phase region except for p = p c = 1/b the specific heat diverges for (/a < 2 as defined 
by a' s and remains finite for the isochore p = p c = 1/b. This demonstrates the exceptional 
character of the critical isochore in this model. 

In the special case that ( = 1 one finds a' s = 2 — 1/a for r < 1 + 2a and a' s = —(3 for 
r > 1 + 2a. Therefore, using ot' s instead of a', the exponent inequalities (|22|) and ( p3|) are 
fulfilled exactly if £ > 1 and r<cr + lorif£ = l and r < 2a + 1. In all other cases (p2|) 
and (|23| ) are fulfilled as inequalities. Moreover, it can be shown that the SMM belongs to the 
universality class of real gases for ( > 1 and r > a + 1. 

The comparison of the above derived formulae for the critical exponents of the SMM for 
C — 1 with those obtained within the FDM (Eqs. 51-56 in |J) shows that these models belong 
to different universality classes (except for the singular case r = 2). 

Furthermore, one has to note that for £ = 1 , a < 1/2 and l-\-a<r<l + 2a the critical 
exponents of the SMM coincide with those of the exactly solved one-dimensional FDM with 



non-zero droplet- volumes |I6| 



For the usual parameterization of the SMM |I] one obtains with ( = 5/4 and a = 2/3 the 
exponents 



at. 
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1 



(28) 



The critical indices of the nuclear liquid-gas PT were determined from the multifragmentation 
of gold nuclei [EH] and found to be close to those ones of real gases. The method used to extract 



the critical exponents (3 and 7' in Ref. [20] was, however, found to have large uncertainties of 



about 25 per cents |21] . Nevertheless, those results allow us to estimate the value of r from the 
experimental values of the critical exponents of real gases taken with large error bars. Using 
Eqs. (17), (19) and (21) one can generalize the exponent relations of Ref. [16] 







and 



r 



2- 



1 + 5 



(29) 



for arbitrary a and (. Then, one obtains with |22| (3 = 0.32 — 0.39 , 7' = 1.3 — 1.4 and 5 = 4 — 5 
the estimate r = 1.799 — 1.846. This demonstrates also that the value of r is rather insensitive 
to the special choice of (3 , 7' and 5. Theoretical values for (3 , 7' and 5 for Ising-like systems 
within the renormalized 4 theory |£| lead to the narrow range r = 1.828 ± 0.001 . 

The critical indices of the SMM were studied numerically in Ref. [I3|. The version V2 
of Ref. |13j corresponds precisely to our model with t = 0, ( = 5/4 and a = 2/3, but their 
results contradict to our analysis. Their results for version V3 of Ref. are in contradiction 
with our proof presented in Ref. fi~IR. There it was shown that for non- vanishing surface 



energy (as in version V3) the critical point does not exist at all. The latter was found in [13] 
for the finite system and the critical indices were analyzed. Such a strange result is, probably, 
due to finite volume effects, although some doubts about the validity of the methods used to 
extract the critical indices (especially, r) remain. 

It is widely believed that the effective value of r defined as t c s = — d\nnk(e)/d\nk attains 
its minimum at the critical point (see references in fl2"l|). This has been shown for the version 
of the FDM with the constraint of sufficiently small surface tension a = for T > T c [25| 
and also can be seen easily for the SMM. Taking the SMM fragment distribution nk{s) = 
g(T)k~ T exp[%k - ^fk a ] ~ k~ T ^ one finds 



r eff = T --k + ^^k° 



r = mm(r eff ) 



(30) 



where the last step follows from the fact that the inequalities a(e) > , v < become equalities 
at the critical point v = a(0) = 0. Therefore, the SMM predicts that the minimal value of T e g 
corresponds to the critical point. 
In the E900 tT 



Au multifragmentation experiment 



the ISiS collaboration measured 
the dependence of T e s upon the excitation energy and found the minimum value min(r e fj) = 1.9 
(Fig. 5 of Ref. [26] ). Also the EOS collaboration [24] performed an analysis of the minimum 
of T eff on Au + C multifragmentation data. The fitted r cff , plotted in Fig. 11. b of Ref. p4 
versus the fragment multiplicity, exhibits a minimum in the range min(r e g) = 1.8 — 1.9 . Both 
results contradict the original FDM ||, but agree with the above estimate of r for real gases 
and for Ising-like systems in general. 



V. Conclusions 



The critical indices of an exactly soluble version of the SMM are derived. The inclusion 
of excluded volume effects generates a principal difference between the SMM and the FDM - 
these models belong to different universality classes. It is found that for the critical exponent 
a! obtained traditionally on the critical isochore Fisher's and Griffiths' exponent inequalities 
are broken. The modification of the definition of a' in the spirit of Ref. (jnj recovers both 
scaling relations. The scaling laws ( p9|) connecting r with /3, 7' and 5 are generalized for 
arbitrary a and (. The ranges r = 1.799 — 1.846 for real gases and r = 1.828 ± 0.001 for Ising- 
like systems are estimated. It is shown that experimental data on nuclear multifragmentation 
agree well with these ranges. Therefore, it is possible that the values for the exponent r < 2 
seen in the data evidence that the critical point of the nuclear liquid-gas phase transition has 
been reached experimentally. 
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